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CREEPING  RAY  ANALYSIS  OF  RESONANCE 
FOR  PROLATE  SPHEROID 


Abstract  --  A  ray  orbit  resonance  condition  is  postulated  and  shown 
to  reproduce  asymptotically  the  resonant  frequencies  of  a  prolate  spheroid. 
The  method  uses  local  creeping  rays  associated  with  a  smooth  convex 
impedance  surface.  As  expected,  the  method  deteriorates  for  high  aspect 
ratio  spheroids  where  the  tip  radii  of  curvatures  are  electrically 
small.  This  work  is  an  alternative  approach  to  Howard's  earlier  geometric 
treatment  and  gives  the  relation  to  creeping  ray  analysis.  The  theory 
provides  impetus  to  the  understanding  of  mode  conversion  and  diffraction 
ray  coupling  near  edges  of  such  scatterers.  Comparison  of  pole  location, 
pole  trajectories  and  layering  with  L.  Marin's  analysis  (1972)^  is  given. 

I.  INTRODUCTION 

There  continues  to  be  strong  interest  in  the  singularity  expansion 
method  (SEM)  and  its  application  to  the  radar  target  identification  problem 
The  singularities  are  the  complex  resonant  frequencies  in  the  u  or  s-plane 
and  exist  in  layers.  Resonance  requires  at  least  one  orbit  of  a  ray  and, 
therefore,  for  the  external  problem,  the  radiation  loss,  associated  with 
the  imaginary  part  of  the  frequency,  is  considerable.  Thus,  the  resonance 
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poles  describe  relatively  late  times  of  a  low  Q  resonator.  This,  is  a 
difficult  regime  for  the  object  identification  problem.  For  this  translates 
to  small  signal  amplitudes  and,  hence,  poor  signal  to  noise  levels. 

In  addition,  the  required  signal  processing  is  non-linear  and  it  happens 
that  small  errors  in  the  signal  can  result  in  large  errors  in  pole  location 
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[1].  To  attempt  to  overcome  these  difficulties,  one  can  take  advantage 
of  the  redundant  information  associated  with  pole  layering.  Thus,  identi¬ 
fication  can  be  based  upon  pole  families,  not  individual  poles.  This 
requires  theoretical  understanding  and  prediction  of- the  ordering  of  the 
higher  order  pole  locations. 

This  is  precisely  the  high  frequency  regime  where  asymptotic  geometrical 
methods  apply.  Indeed,  just  as  classical  mechanics  is  asymptotic  to 
quantum  mechanics  in  the  limit  of  large  quantum  numbers,  so  our  asymptotic 
theory  applies  as  the  discrete  complex  frequency  index  becomes  large. 

Our  geometric  method  also  requires  the  local  radius  of  curvature  to  be 
large  with  respect  to  a  wavelength. 

In  a  previous  paper  by  one  of  us  [2],  the  method  of  Keller  and  Rubmow 
[3]  was  extended  to  the  exterior  resonance  problem.  In  this  paper,  we  make 
explicit  use  of  the  creeping  ray  condition  (see  equation  1),  while  in  the 
previous  paper  [2],  the  equivalent  effect  was  developed  using  an  additional 
explicit  radial  resonance  condition.  This  additional  condition  inforced 
transverse  resonance  between  the  surface  of  the  scattering  object  and  a 
turning  point  caustic. 

Recently,  Heyman  and  Felsen  have  also  shown,  using  a  method  having 
some  similarities  with  ours,  that  local  creeping  rays  describe  exterior 
resonances  on  smooth  convex  objects  [4]. 
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II.  THEORY 


We  make  free  use  of  the  comprehensive  results  of  Bremmer  [5],  Fock 
[6],  Wait  [7],  Logan  [8],  and  more  recently  0.  S.  Jones  [9].  In  these 
references,  it  is  shown  that  the  creeping  wave  mode  equation  for  a  cylinder 
of  radius  p,  and  for  vertically  polarized  waves,  is 

1  v _ 1— —Ji _ !_  =  n  Ml 

(kp)  “"l  (k^p) 

V  ^  I 

* 

In  equation  (1),  is  the  Hankel  function  of  order  v  and  argument 

z  of  type  j  (j  is  1  or  2).  For  a  fixed  frequency  w,  equation  (1)  has  a  discrete 
set  of  roots  v^,  s  =  1,2,3  ....  We  are  using  the  time  dependence  e"^“^.  Here, 
k  =  uj/c  is  the  wave  number  of  free  space  and  k-j  is  the  complex  wave  number 
of  the  spheroid.  Condition  (1)  has  the  physical  interpretation  that  the  sur¬ 
face  reflection  coefficient  is  infinite.  This  allows  a  reflected  wave  with 
no  incident  wave  which  is  a  condition  of  resonance. 

The  propagation  along  the  spheroid  surface  has  angular  dependence 

exp(i  va)  (2) 

where  a  is  an  angular  distance  along  the  surface.  In  our  application,  we 
allow  p  in  equation  (1)  to  be  a  variable.  In  this  ray  treatment,  p  is 
the  local  radius  of  curvature  of  the  spheroid  along  the  direction  of  the  ray 
path.  Hence,  it  is  a  continuous  function  of  position.  Only  the  radius 
of  curvature  in  the  direction  of  propagation  is  significant.  The  transverse 
curvature  is  insignificant.  We  know  that  because  the  local  ground  wave 
modes  for  a  sphere  and  cylinder  of  the  same  radius  and  material  are  equal. 
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(3) 


The  asymptotic  creeping  ray  resonance  condition  is  given  by 

^  V  d  0  =  2Tr(n  +-^)  (3) 

The  term  vd8  is  the  incremental  complex  phase  advance  where  d0  =  —  ,  and 

0 

where  ds  is  the  arc  length  along  spheroid  and  p  is  the  local  radius  of 
curvature.  The  additional  phase  of  tt  corresponds  to  the  phase  shift  advance 
of  ■^/2  going  through  the  focus  at  each  tip  of  the  spheroid.  Details  of 
this  shift  are  given  in  the  appendix. 

There  are  two  geodesic  paths  on  a  spheroid.  First,  there  is  the  path 
from  one  tip  of  the  spheroid  to  the  other;  there  is  a  continuum  of  these. 
They  cross  at  the  tips  which  give  rise  to  the  focus  phase  shift.  Second, 
the  path  going  around  the  equator;  the  resonances  for  this  second  path 
are  easier  to  calculate  because  they  are  the  same  as  the  sphere. 

Returning  to  (3),  along  the  ray  path  implicit  in  the  closed  path 
integral,  we  can  consider  the  geometry  as  two  dimensional.  Along  such 
a  ray  let  the  surface  be  represented  as 


y  =  f(x) 

Then  it  follows  from  elementary  differential  geometry  [10]  that 
de  =  ■-  =  [f"(x)/(l  +  (f'(x))2)]  dx 


(4) 


(5) 


Let  us  apply  this  to  an  elliptic  cross  section  of  a  prol ate  spheroid  for 
principal  orbit  planes  containing  the  axis  of  revolution.  Then 

f(x)  =  ±  (b^  -  c^  x^ 


(6) 


where  c  -  u/ a  i  1  is  the  aspect  ratio  of  the  ellipse  and  where  a  and  b  are 
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the  semi -major  and  semi -mi nor  axes.  Substitution  of  (6)  into  (5)  yields 


de  =  [b^c^/((b^-c^  ci^  x^)  (b^-c^x^))]  dx 


(7) 


where 


The  resonance  condition  (3)  depends  upon  the  propagation  constamt  v.  The 
propagation  constant  v  in  turn  depends  upon  the  local  radius  of  curvature 
j  which  for  the  elliptic  surface  is  given  by 


p  =  (1  +  (f'(x))^) 


3/2 


/f"(x)  =  (b^  -  a^x^) 


3/2 


/(c^b^  ) 


(8) 


In  view  of  the  dependence  of  the  creeping  ray  propagation  constant  v  on 
the  local  parameter  kp,  as  given  by  equation  (1),  it  is  useful  to  make 
a  further  change  of  integration  variable 

=  c^p/b  (9) 


Then  the  resonance  condition  (3)  becomes 


1 


V  dr 


=  (n  +i)  tt/Z 


In  (10)  we  have  accounted  for  the  four  fold  symmetry  of  the  orbit.  Note 
that  the  integrand  in  (10)  is  singular  at  both  end  points  x  =  c  and  1. 
These  singularities  are  removable.  Thus,  we  make  the  linear  change  of 
variable 


(10) 


T 


1  +  c  ,  1  -  c 
2  2 


x 


to  obtain 


(11) 
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2 

1  -  c 


+1 

-1 


X) 


(1  -X^) 
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dx  =  2jr{n  +  i) 


where 


g(u) 


4c  V  (u) _ 

Li(u  +  c)^  (1  +  u)^ 


is  nonsingular. 

Now  we  make  the  further  change  of  variables 


(12) 


(13) 


X  =  sin  Tr/2  s 
Then  (12)  becomes 

f 

g(ly^  +  sin  ( J  s)  ds  =  2iT(n +^)  (14) 

•  -1 

This  form  is  suitable  for  numerical  integration. 

III.  SOME  APPROXIMATIONS 

In  the  high  frequency  limit,  the  creeping  ray  condition  (1)  can  be 
simplified.  The  relevent  approximations  are  discussed  in  detail  in  several 
references  ([5],  [6],  [7],  [8]), 

A  more  modern  mathematical  treatment  is  giv^n  by  Olver  [11].  Olver 
derives  a  uniform  asymptotic  expansion  which  assumes  the  argument  and  order 
are  large  and  approximately  equal.  The  uniform  asymptotic  expansion  is 
given  in  terms  of  Airy  functions.  The  lead  term  in  the  expansion  is  given  by 
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(15) 


hH)  (zv)  =  ^ 


-i  Tr/3  ..  1/4  2/3 

-T/o-  ( - A.  (v  Qai) 

1-Z^  ’ 


where 


2/3  = 


^  ^  (1  -  dt,  w  =  exp  (i  2tt/3) 

1 


(16) 


(17) 


The  phase  factor  w  =  exp(i  2Tr/3)  is  common  in  these  approximations. 

In  fact,  Fock  defined  his  W^(2)  function  as ^ 

(z)  =  -2cj  i  oj^  A.j(ujz) 

This  function  is  outward  radiating  at  infinity.  The  implicit  relation 
(16)  can  be  expanded  as  a  function  of  c.  The  expansion  needed  is  around 
2  =  1  and  hence  5*0.  One  determines 

z(0  =  1  +  3/102"^/^  +  1/700  +  ...  (18) 


Series  (18)  can  be  manipulated  (under  the  same  conditions  required  to 
produce  (8))  to  yield  the  information 


1/3 


v=kp.(^)  t.^V(^)  t^.., 


1/3 


(19) 


2/3 

where  t  =  v  ^  c  and  it  is  assumed  kp  »  1.  In  addition,  series  (18)  can 
also  be  used  to  obtain  the  expansion 


C  =  2^^^  [(1  -z)  +  ^  (1  -z)^  +  ... 


(20) 


2  =  kp/\ 


Also  from  (15)  and  (17)  we  determine 
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h(1) 


(kp) 


2e-i./2  ,-1/3 


4.  1/4 

l-z^)  ' 


where  for  z  =  1 ,  to  first  order,  this  reduces  to 


,(1) 


(kp)  =  2-'-  e 


-  o4/3  ^-i  tt/2  -1/3 


W,(t) 


(21) 


An  explicit  representation  of  the  Fock  function  Wi(t)  is 


W-|(t)  =  -T  exp(st  -  s^/3)  ds  (22) 

7T  *  “(Jj* 

where 

to*  =  exp  (-i  27r/3) 

This  integral  has  two  first  order  saddle  to  points  at  s.  _  =  t/t  which 

5  •  p  • 

coalesce  at  t  =  0.  The  coalescence  corresponds  to  the  position  of  the 
creeping  ray  caustic. 

It  remains  to  determine  an  asymptotic  ray  approximation  tothemodal 
equation  (1). 


IV.  IMPEDANCE  BOUNDARY  CONDITION 

At  high  frequencies,  and  when  the  object  to  which  the  creeping  ray 
attaches  is  convex  and  a  relatively  good  conductor,  it  is  not  necessary 
to  compute  the  fields  on  the  inside  of  the  object.  This  is  possible  because 
under  these  conditions,  the  fields  inside  are  geometrical  and  propagate 
nearly  normal  to  the  surface.  The  ratio  of  tangential  "E  and  H  is  approxi¬ 
mately  the  same  as  for  a  plane  wave.  The  general  formulation  is  known 
as  the  Leontovich  boundary  condition  and  a  good  account  is  by  this  original 
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author  [12].  From  this  discussion  it  follows  that  the  normal  surface  wave 
impedance  in  the  limit  from  inside  is  given  by 


,(1) 


2-k2 


,  horizontal  polarization 


Zi/i^  2_|^2  /k^  ,  vertical  polarization 


(23) 


The  terminology  is  with  respect  to  the  electric  vector.  Thus,  horizontal 
corresponds  to  the  electric  vector  in  the  plane  of  the  surface,  vertical 
corresponds  to  H  lying  in  the  plane  of  the  local  surface.  The  intrinsic 
wave  number  k-j  is  complex,  i.e. 

(<1  =  w/uqS-]  ,  =  e-j  +  i  a/u)  (24) 


Im(ki )  ^  0. 


The  wave  impedance  (23)  is  a  ratio  of  tangential  components  at  the  surface 
of  the  scatterer.  Hence,  the  corresponding  ratio  on  the  exterior  has  the 
same  limiting  value. 


(25) 


On  the  exterior,  it  follows  from  Maxwell's  equations  that 


r  i  w  u  ik/ 


/ 


li 

3n 


( i  w  e  k ) 


,  horizontal 


,  vertical 


(26) 


Substitute  (26)  into  (25).  Let  -k  and  be  the  tangential  E  or  H  field  for 


horizontal  or  vertical  polarization,  where  3/3n  is  the  outward  pointing 
normal  derivative.  Then,  the  Leontovich  boundary  condition  is  explicitly 
given  by 

ii  +  =0  (27) 

3n 

where 


A 


kV'kTT?  /k,' 


,  horizontal  polarization 


,  vertical  polarization 


We  now  relate  this  development  of  the  homogeneous  boundary  condition  (27) 
to  our  creeping  ray  analysis. 

Recall  that  in  equation  (8)  we  have  defined  p  to  be  the  local  radius 
of  curvature  of  the  object  along  the  creeping  ray  path  direction.  Thus, 
the  normal  derivative  3/3n  =  3/3p  .  Then,  using  (21)  and  (27)  yields  the 
explicit  modal  equation 

w.|  (t)  -  q  (t)  =  0  (28) 

where  the  normalized  impedance  parameter  q  is  given  by 


q  =  -i  (v/2)'/3  . 


1/k,  horizontal  polarization 

2 

k/k-j  ,  vertical  polarization 


(29) 


In  the  perfectly  conducting  limit,  the  modal  equation  further  simplifies  to 
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horizontal  polarization 


>  c 


Wi(t')  =  0  vertical  polarization 
I  m 


(30) 


m  =  1 ,  2,  3  . . . 


The  roots  t  and  t'  for  w, (t)  =  0  and  w,  (t)  =0  respectively  are  related 
m  m  I  I 

to  the  roots  of  the  Airy  function  Ai  by  relation  (17).  The  Airy  function 
has  isolated  zeros  along  the  negative  real  axis.  There  are  no  others. 
Thus,  the  roots  defined  by  (30)  all  lie  on  the  ray  i.e. 


m  m 


t;  =  e^' 


(31) 


Abramowitz  and  Stegun  have  listed  the  first  ten  of  them  [13].  They  are: 


m 

m 

^m 

1 

-2.33810741 

-1 .01879297 

2 

-4.08794944 

-3.24819758 

3 

-5.52055983 

-4.82009921 

4 

-6.78670809 

-6.16330736 

5 

-7.9413359 

-7.37217726 

6 

-9.02265085 

-8.48848673 

7 

-10.04017434 

-9.53544905 

8 

-11.00852430 

-10.52766040 

9 

-11.93601556 

-11.47505663 

10 

-12.82877675 

-12,38478837 

Table  1:  Zeros  of  Ai(z)  and  Ai(z) 

(from  Abramowitz  and  Stegun,  p.  478) 

Asymptotically  for  large  m,  the  zeros  a|^  and  aj^  are  given  by 


a^  =  -(3^(4m-l)/8)^''^  +  0  (m""^/^) 

aj^  =  -(37T(4m  -  3)/8)^^^  +  0  (m"^/^)  (32) 
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If  the  surface  is  not  perfectly  conducting,  mode  equation  (28)  must 
be  solved.  However,  Airy's  differential  equation  for  W^(t)  is 

W^'(t)  -  t  W^(t)  =  0  (33) 

This  can  bq  checked,  i.e.  W-j(t)  as  defined  by  equation  (22)  satisfies  (33). 
Thus,  the  roots  t^  as  a  function  of  q  satisfy  the  differential  equation. 


_J _ _ 

“I 

m  ^ 


(34) 


This  is  a  nonlinear  Riccati  equation  which  has  been  solved  in  the  Bremmer 
reference  [6]. 


V.  ASYMPTOTIC  SOLUTION  OF  RESONANT  FREQUENCIES 
Given  the  angular  propagation  constant  v,  (equation  (19)),  the  resonance 
condition  (10)  becomes 


(n  +]■)  J  =  ka  [ 

J  C 


2  . 

T  d  T 


'  (d  -t2)(T^-c2))' 


*  ,2/3  ,^,1/3 


dx 


^  (d  -T^)(t^-C^)) 


.  V3  ,1 


+ _ E  (JL) 

60  'ka^ 


dx 


""  t2((1 -T2)^T2-c2))* 


(35) 


The  integrals  on  the  right  hand  side  of  (35)  are  tabulated  special  functions, 
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The  first  is  k/4  times  the  perimeter  of  an  ellipse  (reference  [14], 
#  223.6c,  p.  49)  i.e. 


rl 

C 


2 

T 


dT 


(d 


(t‘ 


=  E{  /l  -c^  ) 


where 


E(x) 


'Tr/2  2  2  ^ 

(1  -  X  sine)  de 

Jq 


(36) 


is  a  complete  elliptic  integral  of  the  second  kind  (references  [13]  p.  590, 
and  [15]). 

2  2  ^ 

Result  (36)  follows  from  the  change  of  variable  x  *  (1  -t  )/(l  -c  )) 
and  then  the  subsequent  change  of  integration  variable  sine  =  x.  Similar 
changes  of  variables  can  be  made  to  determine 


dx 


^  ((1  --r2)(,2.^2)j^ 


1  +  c 


«Ti|) 


where 


it/2  2  7  “i 

K(x)  =  [  (1  -  x*^  sin'^e)  de 

*0 


(37) 


The  K  function  is  the  complete  elliptic  integral  of  the  third  kind  given 
also  by  Abramowitz  and  Stegun  [13].  The  same  substitutions  used  in  obtaining 
(36),  and  when  an  additional  substitution  x'  =  -  leads  to  the  result 


dx 


- T  ■-^E(/l-C^  ) 

2  .2^,*  c^ 


(38) 
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Mode  equation  (35)  thus  becomes 


K(|^) 


1  +  c 


f/r:?' 


t  X 

m 


(n  )Tr 


4E 


/77? 


+  _!!1_ _  =  0 

I20x  ^ 


(39) 


where  x  =  (ka/2)^'^^, 


Thus,  (39)  is  a  quartic  equation  in  x  with  complex  coefficients. 

3 

In  the  perfectly  conducting  limit,  from  (31)  it  is  seen  that  t  is  real. 
This  suggests  the  change  of  variable. 


X  =  y  t, 


m 


Then  (39)  becomes 

^  1  +  C  I - a- 

e(/i-c2)  tl 


(n+^)  II 


.-2/3 


*  0 


4E(/1 


T2  .6  120  t„y 

c  t„  m  ■' 

m 


(40) 


Notice  the  coefficients  in  (40)  in  the  perfectly  conducting  limit  are  real. 
The  complex  normalized  frequency 


flj^jii  =  ka  ,n,  m-1 ,2,3  . . . 


depends  upon  two  integers,  n  and  m,  which  order  the  roots  in  the  angular 
(e)  and  radial  (p)  coordinates.  In  the  terminology  of  Marin  [16], 

iu)  s 

_  _ _  _  _ nm  _  nm 

"•^nm  ^nm  "  ’  c  ~  c 


(41) 


Q 

where  here  and  what  follows  c  =  2.997925  x  10  m/sec  and  where  is 
the  complex  natural  resonant  frequency  in  radrans/second,  and  s^,^  is  the 
Laplace  transform  pole  location.  In  the  figures  to  follow  it  is  convenient 
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to  plot  the  location  of  the  unitless  frequency 


Y  d 

'nm 


(42) 


where  a  is  the  semi-major  axis  of  the  prolate  spheroid. 

Of  course,  since  the  actual  electromagnetic  time  domain  response 
must  be  real  and  causal,  the  poles  in  the  complex  s  plane  (or  ya  plane) 
must  either  lie  on  the  negative  axis  or  appear  as  complex  conjugate  pairs 
in  the  second  and  third  quadrants.  We,  therefore,  only  plot  results  in 
the  second  quadrant. 


VI.  NUMERICAL  RESULTS 

In  figure  1,  we  reproduce  (from  reference  (16))  a  pole  trajectory  of  the 
first  two  pole  layers  of  the  prolate  spheroid.  Each  curved  line  segment  traces 
the  motion  of  the  poles  in  the  second  quadrant  of  the  normalized  s  plane  as 
the  aspect  ratio  c  »  b/a  varies  from  .1  to  1 .  Figures  2a  through  2d  are  com¬ 
parison  plots  of  Marin's  numerical  solution  of  an  integral  equation  [16]  and 
ours.  As  would  be  expected  of  an  asymptotic  method,  in  the  case  of  a  sphere 
(figure  2a),  the  accuracy  of  our  method  is  good.  When  the  aspect  ratio 
c  =  b/a  departs  from  unity,  the  accuracy  degrades  since  then  at  the  tips  of 
the  prolate  spheroid  the  radius  of  curvature  becomes  electrically  small. 

This  violates  the  asymptotic  assumption.  The  second  layer  plotted  further 
to  the  left  in  figures  2a  through  2d  correspond  to  m  =  2  in  the  ray  equation 
(40).  The  comparison  is  somewhat  misleading  as  tables  2a  through  2d  show. 

The  tables  have  the  same  information  as  the  figures  but,  in  addition, 
include  the  absolute  relate  error  in  percent  entry  in  the  fifth  column. 

The  errors  for  m  =  1  and  m  =  2  pole  layers  are  seen  to  be  comparable. 

Also,  in  these  tables,  it  is  evident  that  the  accuracy  improves  as  the 
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0.2 


I 


Re{ya} 


Fig.  I  The  locus  of  natural  frequencies  when  0.1  <  b/a<I.O.  The  location  of  the 
natural  frequencies  for  b/a  =  0.1, 0.2, 0.5, 1.0  is  indicated  on  the  curves. 
The  arrow  indicates  the  direction  in  which  the  noturol  frequencies  move 
for  increasing  values  of  b/a,  (After  L.  Marin,  1972.) 
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Fig.  2  b  Natural  frequencies  of  a  prolate  spheroid  with  b/a  =  0.5 


Fig.  2c  Natural  frequencies  of  a  prolate  spheroid  with  b/a  =  0.2 
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Table  2a.  Mode  Comparison  for  b/a  = 


b/a  =  0.5 


Integral 
Equation  Method 


Re{ya} 


-  .453 

-  .703 

-  .890 
-1.042 
-1.171 
-1.286 
-1.388 
-1.480 
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-2.879 

-3.408 
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Re{ya} 

Im{ya} 

Error 
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1.232 

24.46 
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2.422 
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Table  2b.  Mode  Comparison  for  b/a  =  0.5 
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b/a  =  0.2 
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Equation  Method 

Ray  Sol 
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Error 

Relya} 

Imlya} 
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(percent) 

-  .336 

1.374 

-  .674 

1.548 

26.88 

-  .516 

2.817 

-  .846 

2.923 

13.78 
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-.-1.075 
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3.98 
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Table  2c.  Mode  Comparison  for  b/a  =»  0.2 


b/a  =  0.1 

Integral 

GTD  Creeping 

Absolute 

Equation  Method 

Ray  Solution 

Relative 

Error 

Relya} 

Imlya) 

ReCya) 

Imlya} 

(percent) 

-  .265 

1.458 

-  .552 

1.704 

25.51 

-  .400 

2.977 

-  .678 

10.61 

-  ,497 

4.510 

-  .772 
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6.051 

-  .848 
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4.39 
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-  ,913 
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-1.022 
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2.46 
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0 
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44.85 
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-2.904 
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32.77 

-4.660 

5.194 

-3.159 

6.202 

25.91 

-4.975 

6.7S2 

-3.382 

7.639 

21.51 

Table  2d.  Mode  Comparison  for  b/a  =0.1 


integer  n  increases.  Thus,  as  also  expected,  our  geometrical  method 
improves  as  either  ka  or  n  increases. 

VII.  CONCLUSION 

We  have  presented  a  geometrical  creeping  ray  formulation  of  the  natural 
oscillations  of  a  prolate  spheroid.  The  general  formulation  has  obvious 
extensions  to  smooth  rotational ly  symmetric  convex  objects.  For  aspect 
ratios  not  less  than  one  half,  the  absolute  percent  error  for  all  but  the 
lower  order  (n  =  1  and  2)  pole  locations  is«less  than  10%.  As  n  increases, 
the  error  becomes  less.  At  the  tip  of  a  prolate  spheroid  of  aspect  ratio 
b/a  =  0.1  |kal  =  10'^  so  clearly  asymptotic  formulas  such  as  (19),  which 
are  fundamental  to  this  method,  are  inappropriate.  In  such  cases,  we 
recommend  that  tip  diffraction  be  implemented.  Indeed,  the  creeping  rays 
loose  their  identity  near  such  tips.  In  other  words,  strong  "mode  mixing" 
occurs.  We  will  consider  this  addition  to  the  theory  in  a  future  paper. 

APPENDIX:  PHASE  ADVANCE  AT  TIP  OF  SPHEROID 

At  a  focus,  it  is  known  that  the  phase  of  the  converging  rays  experience 
a  phase  shift.  Born  and  Wol f  in  their  text  [17]  have  a  section  on  the 
three  dimensional  light  distribution  near  a  focus.  Their  model  is  a  spherical 
wave  emerging  from  a  circular  aperture  and,  hence  conversely,  converging 
on  the  opposite  side  of  the  aperture  toward  an  axial  focal  point.  Although 
geometrical  optics  predicts  a  discontinuous  focus  phase  shift,  finite  wave 
length  calculations  show  it  to  be  continuous.  For  these  optical  treatments 
of  phase  behavior,  the  phase  shift  is  ir  radians.  As  we  shall  now  show 
this  is  not  the  same  as  for  our  nroblem.  At  the  poles  of  a  spheroid,  the 
focus  of  creeping  rays  occur. 
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For  a  specific  model,  we  choose  a  perfectly  conducting  spherical 
surface  of  radius  a.  For  excitation,  we  use  a  vertical  electric  dipole. 
Sommerfeld  gives  a  good  derivation  of  this  model  [18].  After  the  Watson 
transformation,  the  potential  u  is  determined  to  be 


u(r,6)  = 


ik 

4 


+  1 


Sin  v 

-o'^r^2 


7T 


P  (-cos  d) 


kr) 

n  (ka) 


(43) 


In  expression  (43), 

=  (:r/(2z))^  hJ||(z)  (44) 

where  is  the  Hankel  function  of  the  first  kind.  The  denonrinator 

n^(z)  is  given  by 


n  (z)  =  Ch^^^(2)  +  zh^^^  (z)] 

V  3v  V  V 


(45) 


where  prime  denotes  derivative  with  respect  to  z.  In  (43),  ,  m=0,l,2  . 

are  the  roots  to  the  transcendental  equation 


+  zh()_>'(z))  ' 


v=v^,v, .Vj  . 


=  0 
z  =  ka 


(46) 


as  a  function  of  v.  Each  term  in  (43)  is  a  creeping  wave.  For  large  ka, 
Sommerfeld  [17]  shows  that  asymptotically 


V 


m 


ka  [1  +]-(4m+  1)^/^ 


i  ir/3  1 
Mka^  ®  ^ 


(47) 


Since  the  v's  are  large  in  magnitude,  Py(-cose)  can  also  be  approximated 
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by  its  large  argument  asymptotic  expression  e.g.  (see  [17],  p.  149). 

P  (-cose)  =  - ^ - r  cos  (v+4)  (tt-0)  -  Tr/4)  (48) 

(ZTTvsine)^ 

provided  |\;|  >>  1  and  9  not  too  near  0  or  tt. 

Now  the  phase  change  associated  with  the  creeping  ray  passing  through 
the  focus  at  the  poles  of  the  sphere  is  contained  in  expression  (48)  as 
has  been  shown  by  Wait  [  7  ].  We  use  his  argument  and  note  that  result 

4 

(48)  can  be  decomposed  into  two  "rays"  i.e.  we  rewrite  (48)  as 

P  (-cose)  -A  (49) 

V 

where 

-i  (v  +  i)Tr 

A  s  - -  ,•  0  <  9  <  TT 

2(2Trv  sin  9)* 

The  interpretation  is  that  there  are  two  traveling  waves  on  great 
circle  paths.  One  travels  through  an  arc  9,  the  other  through  the  arc 
(2it-9).  Upon  going  through  the  focus  at  9  =  tt,  the  second  ray  path  picks 
up  a  phase  advance  of  tt/2  radians.  Although  not  demonstrated  from  this 
argument,  the  phase  advance  is  not  abrupt  (see  [7])  but  continuous,  and 
Pv(-cose)  near  tt  is  finite  i.e. 

P^(-COSe)  =  Jq((v  +  ^)(ti  -  9))  +  0((Tr-9)^) 


0  =  TT 

(see  [  7  ],  p.  163).  This  completes  our  discussion  of  the  tt/2  phase  advance 
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associated  with  spheroid  focus  phenomena,  and  as  used  in  our  creeping  ray 
resonance  condition(3) . 
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